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Abstract 

A Ginzburg-Landau type equation with nonlocal coupling is derived systematically as a reduced 
form of a universal class of reaction-diffusion systems near the Hopf bifurcation point and in the 
presence of another small parameter. The reaction-diffusion systems to be reduced are such that 
the chemical components constituting local oscillators are non-diffusive or hardly diffusive, so that 
the oscillators are almost uncoupled, while there is an extra diffusive component which introduces 
effective nonlocal coupling over the oscillators. Linear stability analysis of the reduced equation 
about the uniform oscillation is also carried out. This revealed that new types of instability which 
can never arise in the ordinary complex Ginzburg-Landau equation are possible, and their physical 
implication is briefly discussed. 

PACS numbers: 05.45.-a,47.54.+r,82.40.-g 
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I. INTRODUCTION 



Oscillatory react ion- diffusion systems are generally reduced to the complex Ginzburg- 
Landau (CGL) equation by means of the so-called center-manifold reduction when the local 
oscillators are close to their supercritical Hopf bifurcation point This fact led to a now 
widely accepted view that, without resorting to individual systems, one may concentrate on 
CGL if one wishes to gain a qualitative understanding of some universal dyanmical features 
shared commonly by a broad class of oscillatory reaction-diffusion systems. This view in fact 
underlies a vast amount of work in the past (see for review) devoted to CGL. However, 
due to its local nature of coupling, CGL may fail to capture some important aspects of 
the dynamics characteristic to a certain class of oscillatory reaction-diffusion systems. In 
fact, it was argued recently that a situation may arise where the oscillator coupling becomes 
effectively nonlocal, and as a consequence the system exhibits such peculiar dynamics as 
can never be seen in CGL jj, |^ . This suggests that there remains a yet unexplored area 
of reaction-diffusion systems where nonlocal effects on the pattern dynamics must seriously 
be considered. In the present paper, we are concerned with how this area is still accessible 
within the framework of the center-manifold reduction. It will turn out that this can actually 
be achieved by a slight extension of the conventional reduction scheme. 

Effective nonlocality in coupling may become relevant when the reaction-diffusion system 
involves three or more chemical components. Suppose that the system of concern is such 
that the chemical components constituting the local oscillators are free of diffusion, i.e., the 
oscillators remain uncoupled, while the system involves an extra diffusive component which, 
for its diffusive nature, plays the role of a coupling agent. By eliminating mathematically this 
diffusive component, the system becomes a field of nonlocally coupled oscillators possibly 
involving delay also. The main goal of the present paper is to achieve a reduction of such 
react ion- diffusion systems to a universal equation of the Ginzburg-Landau type without 
missing the effects of nonlocality. 

It may seem that the center-manifold idea would not work for our purpose because CGL 
is believed to be quite a general result of the center-manifold reduction. Still there seems 
to be a way out within the same framework if we slightly extend the conventional method 
of reduction. Note first that under usual conditions the disappearance of the effects of 
nonlocality near the bifurcation point comes from that the characteristic wavelength Ip of 
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the pattern becomes far longer than the effective couphng radius. If so, it may happen 
that nonlocahty persists even close to the bifurcation point when the system involves a 
certain parameter whose suitable tuning as the bifurcation point is approached keeps Ip 
comparable with the effective radius of coupling. Similar physical idea lies behind the 
multiple bifurcation theory which aims to capture such complex dynamics as is absent in 
the vicinity of a simple bifurcation point. 

Section II starts with introducing a universal class of oscillatory react ion- diffusion systems 
involving three or more chemical components. After briefly discussing its physical relevance, 
we proceed to its reduction near the Hopf bifurcation. If a certain parameter associated 
with the strength of effective nonlocal coupling is as small as the bifurcation parameter, 
the reduced equation turns out to take the form of a nonlocally coupled complex Ginzburg- 
Landau equation. Analytic formulae for some coefficients of this equation are given, which 
would be of great help in inferring possible ranges of the parameters in the original system 
where nonlocality-dominated pattern dynamics is expected. Regarding the derivation of the 
reduced equation, our primary concern is the case when direct oscillator coupling is absent, 
but the effects of weak diffusive (i.e. direct) coupling will also be considered. In Section III, 
linear stability analysis of our nonlocal CGL is carried out about the uniform oscillation. 
The resulting eigenvalue spectra, especially those of the phase branch, can be qualitatively 
different from those of the standard CGL, whose physical implication is discussed. A short 
summary will be given in the final section. 



II. A UNIVERSAL CLASS OF REACTION-DIFFUSION SYSTEMS AND THEIR 
REDUCTION 

The reaction-diffusion model of our concern was previously proposed by one of the present 



authors 



a 



and is given by the general form 



dtX = fiX)+kgiS), (1) 
rdtS = -S + DV^S + h{X). (2) 

Here the ra-dimensional real vector field / represents a local limit-cycle oscillator with dy- 
namical variable X, so that the first equation with k = represents a field of continuously 
distributed oscillators without mutual coupling. The full system involves an additional chem- 
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ical component with concentration S whose dynamics is governed by the second equation. 
This component simply diffuses and decays at a constant rate, while it is produced locally 
as represented by the term h{X), the production rate depending on the local value of X. 
The dynamics of the local oscillators is influenced in return by the local concentration of S, 
and this effect is represented by the term kg{S). For the sake of convenience, we inserted 
in Eq. (2) parameter r to indicate explicitly the time constant of S, anticipating a limiting 
case in which r is vanishingly small. Note that there is no direct coupling among the local 
oscillators, while their indirect coupling is provided by the diffusive chemical represented by 
S. In later discussions, we will generalize the above model by assuming g and h to depend 
both on X and S , and also by including a small diffusion term i.e. weak direct coupling 
in Eq. (1). Throughout the present paper, the spatial extension of the system is supposed 
sufficiently large. 

Let the spatially uniform steady state of our system be given by {X,S) — (0,0), or 
equivalently, we measure -X" and S always from their equilibrium values. Thus, we clearly 
have h{Q) = 0. It is also convenient to arrange / and g so that the equalities /(O) = g{0) — 
may be satisfied. If we like, S may be eliminated from the system by solving Eq. (2), which 
can be done explicitly because the equation is linear. If the spatial dimension is d, the 
solution of Eq. (2) is given by 

S{r,t) = {271}-'' / dqexp{tqr) / — 

J J -oo T 

exp{-{l + Dq')^-^)hq(t'), (3) 

where hq{t') is the spatial Fourer transform of h{X[r,t')). Equation (1) with S given by 
Eq. (3) constitutes a field of nonlocally coupled oscillators. Note that the nonlocality appears 
in time as well as in space, and the characteristic scales in time and space associated with 
the nonlocahty are given by r and D^/^, respectively. This fact will be used in the discussion 
below. 

Suppose that / involves a parameter fi such that if /i < each local system given by 
Eqs. (1) and (2) with D = has a stable fixed point {X, S) = (0, 0) while this becomes 
oscillatory unstable for /i > 0; namely, at = a pair of complex conjugate eigenvalues of 
the Jacobian associated with each local system about the fixed point cross the imaginary 
axis of the complex plane, while the other eigenvalues all remain in the left half plane. It 
is a well known fact that the center-manifold reduction can be applied to reaction-diffusion 
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systems near the Hopf bifurcation point of the local oscillators This leads generally to 
the complex Ginzburg-Landau (CGL) equation 

dtA = fiaA - (3\A\^A + aV^A, (4) 

where the amplitude A and the parameters a, a and f3 are generally complex. The small 
parameter /i has been retained in Eq. (4) so that A, t and r scale like and 
l/il""*^/^, respectively. 

It is clear that the above reduced form of reaction-diffusion systems remains valid also 
for our particular system given by Eqs. (1) and (2). Since CGL is a diffusively coupled (i.e. 
locally coupled) system, it may seem that effective nonlocality in coupling characteristic to 
our system disappears completely near the bifurcation point. The reason for the disappear- 
ance of nonlocality is clear. This is because the characteristic wavelength Ip of the field X 
becomes longer and longer as the bifurcation point is approached like Ip ~ so that 

the effective coupling radius given by D^^"^ comes to fall well within this scale, which gives 
nothing but the definition of local coupling. In what follows, we will be concerned with a 
special situation in which spatial (and possibly temporal) nonlocality can survive even close 
to the bifurcation point, so that the reduced equation involves nonlocal coupling rather than 
diffusive coupling. The same result was used already in earlier works 0,0] for the particular 
case of vanishingly small r without showing how the reduction can actually be achieved. We 
will develop below the reduction procedure explicitly including the case of finite r. 

Consider our system given by the form of Eqs. (1) and (3) for which Eq. (4) gives the 
right reduced form near /i = provided there is no small parameter other than /i. It is 
clear that the diffusion term in Eq. (4) is the reduced form of the coupling term kg{S) in 
Eq. (1). This implies that |a| = 0(|/i;|). Disappearance of spatial nonlocality is consistent 
with the fact that the characteristic wavelength Ip estimated from the dimensional argument 
for Eq. (4), which confirms Ip = 0(|/u/A;|~^/^), is far larger than the effective coupling radius 
given by D^^"^ provided k remains an ordinary magnitude. This consistency is apparently 
broken if k becomes as small as fi by which Ip becomes independent of fi and hence can 
be comparable with the coupling radius. Putting it differently, spatial nonlocality should 
remain relevant near the bifurcation point provided the strength of coupling between the 
local oscillators and the diffusive component becomes so weak as to satisfy 

A; ~ 0(H). (5) 
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Thus, what we do next is to find a reduced equation vahd near the doubly singular point 
{li,k) — (0,0). Before proceeding to this issue, however, we make a remark on temporal 
nonlocality. Temporal nonlocality, which is characterized by the time constant r, may be 
relevant to the dynamics even near the bifurcation point. However, this effect would never 
appear as a non-Markoffian form of the reduced equation for the amplitude A because the 
time scale associated with the variation of A is much longer than r. Still, r could be 
comparable with the period of the basic oscillation. Then, as we see below, the effect of 
delay may appear in the reduced equation through the change in the space-dependence of 
the coupling function. Specifically, when r is finite, the decay of the coupling function with 
distance becomes oscillatory rather than monotone. 

Suppose first that A; = 0. Then the reduced form of the ordinary equation (1) is given 
by Eq. (4) with a — 0. Although its derivation is routine, we now recapitulate it for the 
purpose of explaining soon later how easily the reduced equation can be generalized when a 
small coupling term is introduced. 

Let the Taylor expansion of /(-X") in terms of X be written as 

f{X) ^ LX + MXX + NXXX + ■■■. (6) 

Regarding the /x-dependence of the coefficients appearing in the above expansion, we need 
to consider it only for the Jacobian L to the first order in n like L = Lq + fJ^Li, higher 
order corrections to L as well as ^u-dependence of M, N etc. are irrelevant to the reduced 
equation to the leading order. Let the pure imaginary eigenvalues at = be ±ia;o, and 
the corresponding right eigenvectors and its complex conjugate are written as U and C7, 
respectively. The corresponding left eigenvectors and its complex conjugate are denoted by 
U* and [/*, respectively. These eigenvectors satisfy U*U = U*U — and U*U — U*IJ — 
1. If is nonvanishing but small, the eigenvalues change to ±iuJo + fiX±, where A+ = A_. 

To the lowest order in yu, the original field X and the complex amplitude A are mutually 
related via 

X{t) ^e'''°^UA{t)+c.c.. (7) 

Thus, in this approximation we have dfA = exp(—iiL!ot)U* {X — LqX). This means fur- 
ther that the right-hand side of Eq. (4) with a = is identical with the reduced form of 
exp{-iLLiot)U*{f{X) - LqX), or 

fxaA - P\A\''A ~ e-^°'U*[f{X) - LqX) 
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= e-''^°*U*{iiLiX + MXX + NX XX + •••). 



(8) 



The standard analysis determines the coefficients a and (3 in terms of some parameters of 
the equations before reduction. It is clear that the linear coefficient a is given by 

a = U*LiU = A+. (9) 

For obtaining (3, the lowest order expression for X given by Eq. (7) is insufficient, and we 
have to use a more precise formula including the next order term: 

X{t) = e^^°*t/A(t) + c.c. 

where 

V+ = V_ = -(Lo - 2iuJo)~^MUU, (11) 
Vo = -2Lo^MUU. (12) 

Using these quantities, /3 is given by 

p = -2U*MUVo - 2U*MUV+ - 3U*NUUU. (13) 

Suppose that the vector field / is modified slightly to f+p. It is clear that the correspond- 
ing reduced equation must also be modified slightly with an additive term exp{—iu!Qt)U*p. 
The specific form of p of our concern is the small coupling term kg{S) in Eq. (1) with 5* 
given by Eq. (3). The original variable X and the reduced one A must now be regarded as 
depending on space as well as on time. Thus, our problem is to find a reduced form of the 
quantity 

kexp{-iu;ot)U*g{S) = kp (14) 

using Eq. (3). Since k is already small, we only need to consider the most dominant contri- 
bution to p. Noting that gr(0) = 0, we may use a linear approximation 

9{S) ^ 9oS 

" (^/^^^"P(^^^)/1V 

exp(-(l + Dq')^-^)hoXq{t'), (15) 
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where — dg/dS\s=o and ho — dh{X)/dX\j^^Q. Thus, using Eq. (7) with X and A 
supposed to depend also on r, we have 



(27r)' 



I dqex.p{iqr) f — 



exp( - (1 + Dq^f—^ - iujoit - t'))Aq{t'), (16) 

where 

rj^{U*g,){hoU). (17) 

Note that Eq. (16) ignores the contribution from the complex conjugate of A{t) which would 
give rise to a rapidly oscillating component of p. This is allowed because such component 
would be averaged out in the equation describing the slow evolution of A{t). Since the 
time-integral in Eq. (16) extends practically over the interval between t — T and t, the slowly 
varying amplitude Aq{t') may safely be replaced with Aq{t). In this way, we obtain 

kp = kr)' J dr'G{r - r')A{r', t), (18) 

where 

Gir) = / ^ge-g^ J + ^ (19) 

^ ^ {27ryJ ^ Dq^ + l + iouor' ^ ' 

and 

V = ■ (20) 

Note that G{t) satisfies the normalization condition 

^ drG{r) = 1. (21) 

Thus the final form of the reduced equation becomes 

dtA = naA - P\A\^A 

+krj' J dr'G{r - r')A{r', t) (22) 

which we call nonlocally coupled complex Ginzburg-Landau equation or simply nonlocal 
CGL. It is clear that the situation of interest is such that k = 0(|/i|) for which the coupling 
term in the reduced equation is balanced in magnitude with the other terms even if the char- 
acteristic wavelength is independent of /i. We assume that the bifurcation is supercritical, 
i.e., the real part of /3 is positive. 
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A few generalizations of the original system (Eqs. (1) and (2)) can be made. Firstly, g 
and h may depend both on X and S. Since the most dominant part of these quantities 
alone is relevant to the reduced equation, one may safely approximate g{X, S) and h{X, S) 
as g{X, S) = g{X, 0) +gi(0, S) and h{X, S) = h{X , 0) + /i(0, S), respectively. The resulting 
new term kg{X, 0) may slightly modify f{X), but the modified f{X) may again be denoted 
by f{X). Furthermore, because S is small, h{0,S) is practically linear in S. Thus, this 
quantity simply modifies the linear decay rate of S which can be normalized by a suitable 
rescaling of time. In this way, the final result of reduction is unchanged except that g{S) 
and h{X) are replaced with g{0, S) and h{X,0), respectively. As the second generalization, 
we may include in Eq. (1) a diffusion term like 

dtX = f{X) + 5 V^X + kg{S) , (23) 

where 5 is a diagonal diffusion matrix with non-negative elements. In parallel with the above 
argument for obtaining a reduced form of the nonlocal coupling term, the reduced form of 
the quantity exp{—iuJot)U*SV'^X will then be added to the right-hand side of Eq. (22). To 
the lowest order approximation, one may apply Eq. (7) for X, by which the above quantity 
becomes SV^A where 5 is a complex number with positive real part, and is given by 

6 = U*6U. (24) 

Thus, Eq. (22) is modified as 

dtA = naA- p\A\^A + SV^A 

+k7]' J dr'G{r - r')A{r', t). (25) 

In the conventional reduction of reaction-diffusion systems, \S\ is assumed to be of ordinary 
magnitude, so that the diffusion term can be balanced with the linear and cubic terms 
in magnitude only if the characteristic wavelength of A scales like However, the 

last property of A contradicts with the particular situation of our concern in which the 
characteristic wavelength remains independent of /i. Therefore, in what follows, we assume 
that \S\ as well as k is of the order of by which all terms on the right-hand side of Eq. (25) 
are balanced with each other, and the coupling nonlocality represented by the last term can 
survive. 
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It is more convenient to write Eq. (25) in the form 



dtA = na'A- p\A\^A + S\/^A 

+ kr]' J dr'G{r-r'){A{r',t)-A{r,t)), 



(26) 



where a' — a + ii~^kr)'. With this form, the couphng term can be approximated by a 
diffusion term when the characteristic wavelength of ^(t",^) is sufficiently longer than the 
coupling radius. Note that a' remains ordinary magnitude because k — 0{ii) by assumption. 

Hereafter, we assume that the system is supercritical or 



An additional remark should be made on the functional form of the coupling function G 
in connection with the time scale r of the diffusive component S. As implied by Eq. (16), 
finite r generally gives rise to memory effects or temporal nonlocality after the variable 
S has been eliminated. However, as was noted before, the reduced equation is free from 
memory effects because the time scale of the slowly varying amplitude A is much longer 
than T. We also noted that because r may be comparable with the period 2t: /ujq of the 
fundamental oscillation, the effect of delay in coupling should be relevant to the reduced 
dynamics. Actually, as is clear from Eq. (19), the effect of finite r appears in the couphng 
function G. For one-dimensional systems, in particular, the coupling function is simply 
expressed as 



and 9 — ujqT. If r is vanishing, G{x) decays exponentially with the decay length being 
given by D^/'^, whereas the decay becomes oscillatory when r is finite. Thus the effect of 
delay in coupling is such that the complex amplitude A{x') in the coupling term is multiplied 
by a factor exp(— ia_|a; — x'\), or equivalently, the phase of A{x\t), denoted by (f){x',t), is 
replaced with (j){x' .f?j — a_|a: — x'\. In the physical language, this means that the phase at 
the spatial point x' experienced by the oscillator at x through the delayed coupling cannot 
be its current value 0(x', t) but should be the value at some time t — tQva. the past, because 
the phase information travels at a finite speed. If this speed is constant, and the oscillation 



Recr' > 0. 



(27) 




(28) 



where 




(29) 
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at x' is nearly regular, then we have (l){x',t — to) equal to (f){x', t) plus something proportional 
to the distance \x — x'\, justifying the above result. 

Coming back to general space dimensions, Eq. (26) with the coupling function given by 
Eq. (19) involves many parameters. However, some of the parameters can be eliminated 
by suitable transformations of some variables. Firstly, the imaginary part of the linear 
coefficient /la' vanishes through the transformation A Aexp[{iiilma')t]. Secondly, one 
may rescale A, t and r in such a way that Re/3, /iRea' and D may all become unity. In this 
way, we may write Eq. (26) as 



dtA = A-{l + ic2)\A\^A+{Si+iS2)V^A 
+K{l+ici) Jdr'G{r - r'){A{r', t) - A{r, t)), (30) 

where the coupling function G is defined as an integral form given by Eq. (19), or 

'^<''> = (2^ 

The reduced equation now involves six independent parameters ci, C2, K, 9, Si and S2 
all of which are independent of smallness parameters /i, k and S. Note that the coupling 
coefficient in Eq. (30) is related to some original parameters through 

K{1 + ^cl) = = ^ ^"^'^ , (33) 

or 

Reo" , , 

K = l- —— 34 

Rea' ^ ^ 

so that, by combining the last equation with the inequality (27) and the original assumption 
Recr > 0, we have a restrictive condition 

K <1. (35) 



III. EIGENVALUE SPECTRUM ABOUT THE UNIFORM OSCILLATION 

It is clear that Eq. (30) admits a family of plane wave solutions Ak{r, t) — Rexp[i{kr — 
fit)] the stability of which is extremely important to the understanding of the pattern dy- 
namics of our system. In the present paper we will focus on the stability of the uniform 
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oscillation Ao{t) for which 



R ^ 1, (36) 

<^o = -C2- (37) 

We now put 

A{r,t)^{l + A{r,t))Ao{r,t), (38) 

and linearize Eq. (30) in A(r,t). The linearized equation can be solved in terms of the 
Fourier components of A{r,t), denoted by Ag(i). Assuming its time-dependence in the 
form A.q{t) oc exp(Ai), we find the eigenvalue equation 

+ u{q)X + v{q) ^0, (39) 

where 

u{q) = -2Re7(g), (40) 
v{q) = hiq)\'-m\', (41) 

with 

7(g) = -{I+IC2)- i6i + i6.2)q^ 

+K{l + ici){G,-Go). (42) 

In what follows, we shall first concentrate on the case without diffusive coupling, i.e. 61 = 
S2 — 0; the effects of nonvanishing but small Si will be touched upon later. 

Let the solutions of Eq. (39) be denoted by A+ and A_ with ReA+ > ReA_. The uniform 
oscillation is stable if and only if ReA+ < holds for all q. This holds only when u{q) > 
and v{q) > for all q. Equivalently, if one of these inequalities becomes violated for a certain 
q, then the uniform oscillation loses stability. This implies that the types of instability could 
be classified in terms of the signs of u{q) and v{q). Simple calculation shows that u{q) and 
v{q) can be expressed in the following form where we use the notation Q = q^. 

n{q) = aQMQ), (43) 
v{q) = aQ)UQ)Q- (44) 
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Here CoiQ) and Co{Q) are non-negative functions of Q, and 





= a2(5^ + aiQ + ao, 


(45) 


CiQ) 


= hiQ + 60, 


(46) 


n f\ 




v^' ) 


Oi 


= 2 + ir(l + ci^), 


(48) 


02 




(49) 


&0 


= 2ir(l + CiC2 + (Ci-C2)e), 


(50) 


&1 


= ir(2(l + ciC2) + ir(l + c?)). 


(51) 



Figure 1 shows the ai-a2 plane divided into three domains (labeled by A, B and C) corre- 
sponding to qualitatively different forms of u{q). Note that the 01-02 plane covers all possible 
types of u{q) because oq is positive definite. Similar picture for v{q) in the h^-hi plane is 
displayed in Fig. 2 where the whole space is divided into four domains (labeled by A' to 
D'). Since the parameters oi, 02, 60 and 61 as functions of four independent parameters K, 
6, c\ and C2 can be changed independently, every combination between two members, one 
from the group (A,B,C) and the other from (A',B',C',D'), is possible. Uniform oscillation is 
unstable for all these combinations except for (A,A'). Loss of stability occurs as we move 
across the line 02 = or — 4aoa2 = in Fig. 1, for which the instability is oscillatory, 
or otherwise 60 = or 61 = in Fig. 2, for which the instability is non-oscillatory. The 
critical wavenumber for which ReA < becomes first violated equals zero on the line 60 = 0, 
finite on q\ — 4aoci2 = and infinite on 02 = and 61 = 0. The last type of instability, i.e. 
the instability which starts at infinite g, was called short-wavelength bifurcation by Heagy 
et alfl. 

If \K\ is not too large, u is non-negative, so that the instability is only through v becoming 
negative and hence it is non-oscillatory. The corresponding critical lines = and bi = 
in Fig. 2 can now be translated into critical relations between K and Kci under fixed C2 
and 6. In this way, we have four types of eigenvalue spectrum shown schematically in Fig. 3 
where lines Li and L2 correspond to 60 = and bi = 0, respectively. The figure shows 
the spectra only for phase-like fluctuations, because in each case the amplitude branch 
which is separated from the phase branch remains negative for all q so that the amplitude- 
like fluctuations are irrelevant to stability. Separation between the phase and amplitude 
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branches also implies that the eigenvalues associated with these branches are all real. In 
contrast, if these branches merge, then the eigenvalues of the two branches for given q 
would form a complex conjugate pair. Note that in each type of spectrum shown in the 
figure the eigenvalue saturates to a constant as the wavenumber tends to infinity, which is 
characteristic to nonlonally coupled systems. The characteristic wavenumber about which 
the eigenvalue starts to saturate equals the inverse of the coupling radius. Physically, this 
reflects the fact that the dynamics of fluctuations whose wavelengths are much shorter than 
the coupling radius is practically the same as the oscillators' individual dynamics, and hence 
it is insensitive to the wavenumber. Note also that the type U2 spectrum is possible only 
when 9 is nonvanishing or, equivalently, when the nonlocal coupling before reduction involves 
delay. 

For larger \K\, the phase and amplitude branches can merge, so that oscillatory instability 
becomes possible. In particular, the short-wavelength type instability i.e. the instability 
initiated by fluctuations with infinite wavenumber can occur across the line 02 = or 
K — —1. This type actually appears in the stabihty diagram in Fig. 4 which is a global 
extension of the diagram in Fig. 3. The values of the other parameters are the same as in 
Fig. 3, i.e., C2 = 1.0 and 9 = 6.0. Oscillatory instability initiated by fiuctuations with finite 
wavenumber, which occurs when crossing the line af — 4aoa2 = in Fig. 1, does not appear 
in Fig. 4, but may appear when different values of C2 and 9 are chosen. For instance, the 
stabihty diagram for the case of C2 = — 2.5 and ^ = 1.5 is displayed in Fig. 5 where the fine 
L4 gives the boundary associated with the last type of instability. 

Finally, we comment on the effects of weak diffusion represented by the term {6i+i52)V^A 
in Eq. (30). For simplicity, we assume that 62 is vanishing while Si is a small positive. It 
is clear that this diffusion term simply gives rise to an additional stabilizing term — to 
each of the eigenvalues A± as a function of q. Thus, the fluctuations with sufficiently short 
wavelength are always decaying. A particularly interesting result from this fact appears in 
type U2 eigenvalue spectrum in Fig. 3. Depending on the value of 5i, the dispersion curve 
is deformed like some curves shown in Fig. 6. Even the dispersion curves of types Ui and 
U12 could be deformed into a similar form if some parameter vaules are chosen suitably. It 
is clear that there is a critical value of Si at which fluctuations with a certain wavenumber 
qc become unstable. Due to the presence of long-wavelength fluctuations which are almost 
neutral in stability, the unstable growth of the mode with qc will generally be unable to 
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lead to a Turing type periodic standing pattern. Instead, from the outset, a group of modes 



with wavenumbers close to qc will couple nonlinearly with another group of modes of &. 



vanishing q with almost neutral stability, leading to a peculiar spatio-temporal chaos 
The simple evolution equation 







most 



121. 



d,u = -dl[e - {I + dlf]u - {d^u)\ (52) 

called the Nikolaevskii equation j^, gives quahtatively the same eigenvalue spectrum as those 
in Fig. 6. Actually, it was argued previously that the Turing pattern in this system existing 
for small positive e is always unstable, and as a result the system immediately becomes tur- 
bulent characterized by the coexistence of turbulent fluctuations with vastly different length 
scales. Possibility of this type of turbulence in react ion- diffusion systems was discussed by 
Fujisaka and Yamada^^. In electro-convective systems, similar type of complex behavior 
was discovered by Kai et al.j^ which they called soft-mode turbulence. 



IV. SUMMARY 



If spatially distributed limit cycle oscillators are diffusively coupled not directly but via 
a certain diffusive component, the oscillators may be viewed as coupled nonlocally after a 
mathematical elimination of this diffusive variable. We showed in the present paper that this 
actually occurs in reaction-diffusion systems and that the nonlocality of this kind persists 
even close to the Hopf bifurcation point provided the coupling between each local oscillator 
and the diffusive component is sufficiently weak. Specifically, under this condition, the 
systems is reduced to a comlex Ginzburg-Landau type equation with nonlocal rather than 
diffusive coupling. Temporal nonlocality or memory effects which generally exist before 
reduction does not appear explicitly in the reduced equation, still they may generally affect 
the functional form of the coupling function. Our results were generalized so as to include 
the effects of direct but weak diffusive coupling among the oscillators. 

Linear stability analysis of the uniformly oscillating state of the reduced equation was 
also carried out. Some new types of eigenvalues spectrum were found to arise, and their 
physical relevance was suggested. 

How the solution behaves in the nonlinear regime is known only partially. For instance, 
type Ui2 spectrum in Fig. 3 was found to lead to turbulence with multi-affinityj^. It was 
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also found that even the normal spectrum (type S) can give rise to peculiar spiral waves 



without phase singularity when the nonlocal coupling becomes weak 



ll|. Different aspects 



of the nonlocal CGL including the afore-mentioned type of turbulence associated with the 



dispersion curve in fig. 6 will be developed in forthcoming papers 



121. 
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Figure Captions 

Fig.l The sign of u{q) changes with q in three different ways. The figure shows how the 
corresponding domains A, B and C appear in the 01-02 plane. u{q) is positive for all 
q in A, negative above some critical q in B, and negative in a finite interval of q in C. 
The line separating domains A and C is given by a parabola 02 = al/AaQ. 

Fig. 2 The sign of v{q) changes with q in four different ways. The figure shows how the 
corresponding domains A', B', C and D' appear in the 60-^1 plane. v{q) is positive for 
all q in A', negative above some critical q in B', negative for all q in C ', and negative 
below some critical g in D'. 

Fig. 3 Four types of eigenvalue spectrum for the phase-like fluctuations about the uniform 
oscillation which is stable only for type S. The figure shows how they appear in the 
K-Kci plane in the vicinity of the origin K = Q. The other parameters are fixed as 
C2 — 1.0 and 6* = 6.0. In each of the four cases, the eigenvalues of the amplitude-like 
fluctuations, which is not shown in the figure, form a branch completely separated 
from the phase branch and remain negative for all q. 

Fig.4 Figure 3 is extended to the region of larger K. A new type of eigenvalue spectrum 
appears for large negative K. This type, labeled by U3 is similar to U2 in Fig. 3 except 
that the amplitude branch merges with the phase branch above a certain wavenumber 
which is still below the wavenumber at which ReA changes sign. Thus, the eigenval- 
ues associated with the unstable fiuctuations are complex rather than real, i.e., the 
instability is oscillatory in nature. 

Fig. 5 Similar to Fig. 4 but for different values of C2 and 9, i.e., C2 = —2.5 and 9 = 1.5. New 

type of eigenvalue spectrum, denoted by U4, appears which differs from U3 in Fig. 4 
only in that the eigenvalues saturate to a negative value as q goes to infinity. 

Fig. 6 Eigenvalue spectra close to a finite-wavenumber instability, which are obtained as a 
modification of U2 in Fig. 3 by assuming nonvanishing 5i in the reduced equation (30). 
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